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Abstract 

Classic and recent results for the critical behaviour of ideal Bose gas at con- 
stant volume and constant pressure and for various spatial dimensionalities d > 
are reviewed. New results about the critical properties in a close vicinity of the 
A— point are presented. 
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1 Introduction 

The thermodynamic potential and the correlation functions of ideal Bose gas (IBG) PjJ 
12] were studied in a number of preceding papers [13 IU 03 EJ IE1 EH UD] ; see also the 
review |TT]. In this paper we shall present a brief review of known results together with 
an unified treatment of critical properties of IBG based on the derivation of effective 
free energy of Landau- Ginzburg type. Our consideration makes possible to outline a 
general picture of the critical behaviour of IBG as well as to present some new results 
about the effects of thermodynamic constraints: the conditions of constant volume 
(V = const) and constant pressure (P = const) for certain spatial dimensionalities 
(d > 0). We should emphasize that the thermodynamic behaviour of IBG near the 
A— point corresponds either to a phase transition of first order or, in cases which we 
shall enumerate, resembles a continuous phase transition that is not equivalent to the 
standard second order phase transition described by tp A — model (see, e.g., Ref. |2J). Our 
results could be used in the interpretation of thermodynamic behaviour of extremely 
dilute Bose-Einstein condensates (BEC), where the interparticle interaction can be 
neglected. 
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2 Thermodynamic potential 



We consider noninteracting and homogeneous gas of spinless bosons, described by the 
Hamiltonian IH EJ 



k 

where r = — fi > is the modulus of the chemical potential fi, k = {kj\j = 1, ...,d} 
is the (i-dimensional wave vector, and at are second quantization operators, and 
iV is the total number of bosons. The external field h generates the spontaneous 
symmetry breaking ((So) > 0) below the critical temperature T c . The energy spectrum 
€k = ck a includes two cases: (i) a — 2, corresponding to real boson particles and boson 
excitations produced by short-range interparticle (fermion and spin) interactions, and 
(ii) < a < 2, which stands for boson excitations (composite bosons) produced by 
long-range interparticle interactions (see, e.g., Refs. [21 E])- The parameter c is equal 
to (h 2 /2m), where m is the effective mass of the bosons. 

The grand canonical thermodynamic potential of IBG, 

Q(T, r, h) = -/T 1 In Tr exp , (2) 

can be exactly calculated [3 lU E] • This quantity can be written in general form 



1 hh* 

n(T, r, h) = -p-WX^A^d, cr)g d/a+1 ((3r) - - — , (3) 

JM r 

where Xt = (2irh 2 /raksT) 1 /' 7 is the thermal wavelength, 



is the integral Bose function (see, e.g., j2J) and 

2 l- d+ 2 d /, T{d/a) 

A ^ a >- ff7r rf(i/2-i/ ff )r(d/2) [ ) 

is a constant with the property A(d, 2) = 1. Note, that the grand potential (3) obeys 
differential relation dQ(T, r, h) = —SdT + Ndr — ^>dh* — ty*dh, where S is the entropy 
and (\& ~ (do)) is the order parameter of BEC. 

The free energy Q(T,r, \1/), in which the natural thermodynamic variable is the order 
parameter can be obtained from Q(T, r, h) through the standard Legendre transfor- 
mation 



fi(T, r, *) = n(T, r, h) + h** + h**\ m , 
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(6) 



where 



on 

dh* 



h 



■iff* = 



on 

dh 



h* 
iVr' 



(7) 



Using Eqs. (3)- (7) we obtain 



fi(T, r, *) = -/T VA^A(d, (r)g d/a+1 (Pr) + Nr^ 2 . (8) 

We can restrict our calculations without a loss of generality to real values of h and 
The susceptibility is given by xt = d^ /dh = 1/Nr ~ t -7 , where 7 is the susceptibility 
critical exponent and t = (T — T c )/T c (see, e.g., [2]). 

For small r and energies (e^ <C /cbT) the correlation function x{k) = fiN~ l (a^(i-g) 
takes the form A^ _1 (c/c <J +r) _1 . This gives the correlation length £ = (c/r) 1 / " ~ t -1 ' and 
~ at r = 0. From ~ ^ _<J we obtain that the Fisher exponent n, defined by 
x{k) ~ k~ 2+v , is equal to (2— a) for all dimensional ranges and constraints (see Table 1). 
In order to obtain the correlation length exponent v and the exponent 7 of susceptibility 
Xt, we need to derive the function r(t) for a general spatial dimensionality d for cases 
of constant volume and constant pressure. The nullifying of chemical potential defines 
the critical temperature T c of BEC for both cases (V = const, P = const). 



3 Constant volume 



For the thermodynamic condition of constant volume we need the free energy 
f(T, N/V,ty), in which the natural variable is the volume V (or the number density 
N/V). To obtain it we expand the integral Bose function g u (y) in powers of fir and 
express the chemical potential as a function of the volume (or density). We expand 
the function g v in the potential Q to the same order in fir. This allows to make the 
next Legendre transformation 



f(T, N/V, *) = uj(T, r, *) - ^r\ r(N/v) , (9) 

where 



y = ^ = ^T d A(d,a)g d/a Wr) + y^\ 



(10) 
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Figure 1: Effective free energy for d < a (d = 1, a = 2), V = const: line 1 (T = 0.25), line 2 
(T = 0.15), line 3 (T = 0.02). The case d = a is qualitatively the same. 

f = F/V and uj = Q/V. Subtracting the free energy of the disordered phase fd is 
we obtain the following effective free energies corresponding to several dimensionality 
ranges: 



ma/{a-d) 

f(T,N/V,*)~- — 

y/(a-d) d < a 

( 1 - fr 2 N \ 
~ 1 exp — j d = a 

J \ d/{d-a) 

-m 2 1 I o < d < 2a 

a 

(fr 2 + 2t) 2 -. - ' — d = 2a 

V ; ln(^ 2 + 2t)- 1 



, 2 

^ ■ ■■' 2 



^ Z + t) d>2a. (li; 

The case d < a is depicted in Fig. 1 for d — 1 and o = 2. We take the volume 
per particle v = V/N to be unity and m — K — ks — 1- The number coefficient in 
Eq. (11) is equal to 1/2 and T c = 0. When the temperature tends to zero the curve 
f{^) acquires a flat part, indicating that the order parameter \l/ occurs with a jump at 
T c = 0. In Fig. 2 we show the d > 2a case, when d = 5 and o = 2. Now the numerical 
factor in Eq. (11) is 8.96 and the normalization ^ = ^Ja/dfy has been used. The 
function /(^), shown in Figs. 1 and 2, is quite different from the respective potential 
/( 1 J r ) known for the usual cases of IBG (d — 3, a = 2) jl] and the y? 4 — model j2j. 

Using the free energy (11) we calculate the specific heat at constant volume Cy = 
—T(d 2 f/dT 2 ) and the equation of state h = df/dty. Bearing in mind the definitions 
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order parameter 



Figure 2: Effective free energy for d > 2a (d = 5, a = 2), V = const: line 1 (t = 0.5), line 2 (t = 0), 
line 3 (t = —0.4). The case d = 2a is qualitatively the same. 

(see, e.g., 0), C v ~ * ~ (A ~ 0) and * ~ ^ (h = 0) we obtain the critical 
exponents a, 6 and (3. The values of the critical exponents are given in Table 1 . The 
behaviour of V, S, Cy and \t at T c is given in Table 2. 



4 Constant pressure 

When the system is under constant pressure we should use the free energy per particle 
fi, in which the natural variable is the pressure P . To obtain the chemical potential 
we solve the equation for the pressure P = —df/dv with respect to the volume per 
particle v to the lowest order in t. Thus we make the following Legendre transformation 

/i(T, P, *) = f(T, N/V, *) + Pv\ v(P) , (12) 

where [L— &/N and / = F/N = fv. 

For the various dimensional ranges we obtain 

1 In Table 1 we distinguish logarithmic factors by "exponents" with suffix I. Thus fit) is said to 
have exponent xi or x , respectively, if f(t) ~ ^/Int^ 1 or f(t) ~ i^lni -1 0. 
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/i(T,P,^)~-(l-^ 2 )W d (l + at) d<a 



~ -(1 - * 2 )t(l + at d/a - 1 ) a<d<2a 
~ -(1 - ^ 2 )t + (1 - ^ 2 )t 2 (l + at) \nr l d = 2a 

~ -(l-^ 2 )t-(l-<I> 2 )t 2 {l + at d/a - 2 + bt) d>2a. (13) 



In Fig. 3 we depict the case a < d < 2a for d = 3 and a = 2. We set again 
m = h = ks = 1 and take the magnitude of the pressure P so as to give T c — 1. 
This gives the factor 1.28 for the coefficient of proportionality in Eq. (13). This picture 
should be compared with those shown in Figs. 1 and 2 (see, also, [21 II])- 

We derive from the free energy fi the specific heat at constant pressure Cp = 
—T(d 2 fi/dT 2 ) and the exponent a. Moreover the chemical potential (13) can be used 
to calculate the correlation length and susceptibility exponents v and 7. Since the or- 
der parameter has a jump, as is seen from Eq. (13) (/x is zero below T c ), the exponents 
j3 and 5 cannot be defined. The values of the critical exponents and the behaviour of 
V, S, Cp and xt at T c are given in Tables 1 and 2. 
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5 Summary 



In this paper we derived and analyzed the IBG free energies of Landau-Ginzburg type 
for the thermodynamic constraints of constant volume and constant pressure at dif- 
ferent spatial dimensionalities d > 0. The values of the critical exponents for various 
spatial dimensionalities are given in Table 1. The critical exponents for d = (a, 2a) and 
P = const are new results of our investigation. The derived free energies enable us to 
present the behaviour of thermodynamic functions at the critical point. It is shown in 
Table 2. Besides, our results make possible to estimate the order of the phase transition 
to BEC. When d > a and P = const this phase transition is of first order. But for 
lower dimensionalities (d < a) and in both cases (P = const, V = const) , the phase 
transition to BEC resembles a first order phase transition only by the order parameter 
jump. On the other hand, for d > a and V = const the phase transition to BEC 
is a continuous phase transition with properties, which are quite different from the 
properties of the standard second order phase transition described by y? 4 — model. Our 
results can be used in the interpretation of experimental results for BEC in extremely 
dilute boson gases at constant pressure or volume. 



Table 1. Critical exponents of the ideal Bose gas. 



VOLUME=COnst 


a 




p 


7 


5 


V 




d<a (T c = 0) 


-d/a 




* 


a 
c—d 


* 


l 

a— d 


2- a 


d = a (T c = 0) 


-1 




* 


OO 


* 


OO 


2-a 


a < d < 2a (T c > 0) 





(2-2(7 
d— a 


1/2 


(7 

d—a 


d+a 
d—a 


1 

d—a 


2-a 


d = 2a (T c > 0) 







1/2 


h 


3* 


av = l t 


2-a 


d > 2a (T c > 0) 





^2a— d 


1/2 


1 


3 


1/a 


2-a 


PRESSURE=COnst 


a 




P 


7 


5 


V 


V 


d < a (T c > 0) 


2 - a/d 


-l b 


* 


a/d 


* 


l/d 


2-a 


d = a (T c > 0) 


0i 




* 


h 


* 


av = li 


2-a 


a < d < 2a (T c > 0) 


2 -d/a 




* 


1 


* 


1/a 


2-a 


d = 2a (T c > 0) 







* 


1 


* 


1/a 


2-a 


d > 2a (T c > 0) 





(^rVi) a 


* 


1 


* 


1/a 


2-a 



*The order parameter has a jump (from to 1) at T c . 
a The first entry is for d < 3cr, the second is for d > 3a. 
b This entry is only for d = a/2. 
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Table 2. Behaviour of *, V, S, C v , C P , xt at T c . 



VOLUME=COnst 


order par. 


volume 


entropy 


heat cap. 


suscept. 


d<a (T c = 0) 


jump 


const. 


cont. 


cont. 


diver g. 


d = a (T c = 0) 


jump 


const. 


cont. 


cont. 


diver g. 


a < d < 2a (T c > 0) 


cont. 


const. 


cusp 


cusp 


diver g. 


d = 2a (T c > 0) 


cont. 


const. 


cusp 


cusp 


diver g. 


d > 2a (T c > 0) 


cont. 


const. 


cusp 


cusp 


diver g. 


PRESSURE=COnst 


order par. 


volume 


entropy 


heat cap. 


suscept. 


d < a (T c > 0) 


jump 


cont. 


cont. 


cont/jump/div. a 


diver g. 


d = a (T c > 0) 


jump 


cont. 


cont. 


cont. 


diver g. 


a <d<2a (T c > 0) 


jump 


jump 


jump 


divergent 


diver g. 


d = 2cr (T c > 0) 


jump 


jump 


jump 


In divergent 


diverg. 


d > 2cj (T c > 0) 


jump 


jump 


jump 


jump 


diver g. 



"Respectively for d < a/2, d = a/2, d > a/2. 
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